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SPACE- CHARGE -INDUCED  BREAKDOWN  IN  DIELECTRICS 

Massachusetts  Institute  of  Technology 
Laboratory  for  Electromagnetic  and  Electronic  Systems 
High  Voltage  Research 


1.  PROJECT  OVERVIEW 

This  repon  describes  the  results  of  a  15  month  study,  sponsored  by  AFOSR,  Directorate  of 
Physical  and  Geophysical  Sciences,  concerning  dielectrics  and  more  particularly  the  influence  of 
accumulated  space  charges  on  dielectric  breakdown.  The  work  includes  experimental  tests  and 
theoretical  modelling  of  space  charge  effects.  The  tests  employed  energetic  electron  beams  as  a  means 
to  inject  controlled  quanrides  of  charge  into  test  specimens.  The  use  of  electron  beams  is  relevant  to 
cenain  pracdcal  applicadons  such  as  radiadon  effects  and  space  craft  charging,  as  well  as  being  an 
important  diagnosdc  tool  for  dielectric  materials  and  for  e-beam  transport  smdies. 

The  invesdgadons  were  organized  into  5  smdy  tasks:  discharge  extent,  charge  transfer, 
channeling  locadon,  flux  rate  effects,  and  a  propagadon  model.  The  results  firom  these  tasks  cover 
many  important  issues,  which  for  purposes  of  this  report  have  been  grouped  into  the  broad  categories 
of  measured  characterisdcs  and  models.  Specifically,  the  repon  includes,  after  the  introducdon  (2), 
secdons  on  the  characterisdcs  of  induced  breakdown  (3),  modelling  of  the  charging  process  (4), 
modelling  for  breakdown  within  a  charged  volume  (5),  new  diagnosdc  methods  (6)  and  conclusions 
(7). 

2.  INTRODUCTION 
2.1  Space  Charge  Effects 

Dielectrics  by  defiidtion  impede  the  flow  of  electric  charges  and  hence  can  act  to  store 
electric  charges  on  their  surface  or  in  their  volume  and  thereby  cause  major  changes  in  electric  field 
distribudons.  The  residence  time  for  such  stored  charge  varies  from  material  to  material,  with  the  type 
of  carrier,  and  according  to  the  type  of  transport  process.  When  the  accumulated  space  charge  is  large, 
it  can  produce  electric  fields  which  can  exceed  those  which  the  dielectric  can  sustain  and  breakdown 
can  occur. 


Situations  where  space  charges  are  often  important  include  a  range  of  di^erent  applicadons 
and  can  be  separated  into  cases  with  and  without  e-beams.  When  e-beams  are  not  present,  space 
charges  can  still  accumulate  but  by  other  charge  transport  processes. 
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Example  applications  include: 

With  e-beams: 

-  basic  beam  interaction  physics 
-dosimetry 

•  ladiadon  damage 

-  multichannel  switching 

-  accelerator  technology 
Without  e-bcams: 

-  electrical  treeing 

-  aging,  degradation 

-  reliability 

-  surface  fiashover 

-  solid  dielectric  interfaces 

Models  to  describe  the  charge  accumulation  stage  are  usually  based  upon  the  continuity 
equation  which  is  coupled  to  Poisson's  equation  for  the  electric  field.  A  driving  term  which  accounts 
for  charge  injection,  such  as  from  an  e-beam  is  needed  to  activate  the  process.  These  equations  are: 

7  .  j  +  l£.  ,  0 

3t 

J  *  oE  +  ppE  +  Jq 
V  .  E  »  p/e  ,  /  E  .  dl  »  -V 

where  J  is  the  current  density,  p  the  net  charge  density,  o  the  material  conductivity,  E  the  local 
electric  field,  the  carrier  mobility,  Jq  the  current  density  due  to  external  generation  sources,  and  c 
the  dielectric  permittivity. 

2.2  Basic  Breakdown  Processes 

When  space  charge  generated  fields  are  high  and  breakdown  occurs  it  is  a  sudden 
irreversible  event  with  channels  which  rapidly  develop  and  propagate  through  the  dielectric.  Space 
charge  induced  breakdown  has  some  characteristics  which  are  in  common  to  those  of  usual  applied 
voltage  breakdown  and  in  addition  some  characteristics  which  are  unique.  As  developed  in  this  smdy 
space  charge  effects  were  isolated  by  nai  using  an  electrode  pair  to  apply  an  electric  field  to  the 
dielectric.  Without  an  externally  applied  field,  only  the  accumulated  space  charges  held  by  the 
dielectric  produced  the  electric  fields  and  hence  space  charge  effects  could  be  clearly  isolated  for  smdy. 
The  conditions  associated  with  charge  accumulation,  channel  initiation  and  channel  propagation  were 
of  interest  in  this  work. 
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Models  for  the  breakdown  process  arc  less  developed  than  those  for  the  charge 
accumularion  stage.  However,  by  use  of  controlled  distributions  of  implanted  charges  as  employed  in 
this  work,  it  has  been  possible  to  determine  major  parameters  of  importance  to  the  breakdown  process. 

2.3  Dielectric  Materials 

The  role  of  the  dielecmc  material  on  space  charge  induced  breakdown  involves  three  main 

processes: 

-  accumuladon 

-  breakdown  incepricn 

-  breakdown  extent,  treeing. 


Firstly,  the  material  influences  the  net  charge  accumulation  by  its  inherent  charge  transport  and 
conduction  processes.  Second,  once  charge  has  accumulated,  basic  material  insulation  strength  and 
local  defects  influence  the  initiation  of  breakdown  channels.  And  thirdly,  channel  propagation  is 
influenced  by  chemical  and  plasma  details  including,  power  loss  rates,  plasma  waves,  and 
conductivities  in  the  channel.  In  these  investigations  many  of  the  tests  were  made  with 
polymethylmethacrylate  (PMMA)  since  it  readily  accumulated  charges  and  is  transparent  for  optical 
analysis.  Other  materials  included  polycarbonate,  epoxy  and  Mylar. 

2.4  Experimental  Methods 

These  studies  have  employed  an  unusual  procedure  for  insulation  studies,  the  direct 
deposition,  ie.  implantation,  of  charges  within  the  test  dielectrics.  Macroscopic  volumes  of 
space-charge  were  implanted  within  the  solid  dielectric  samples  by  means  of  an  energetic  electron 
beam.  To  do  this  a  Van  de  Graaff  accelerator  with  controlled  beam  energies  up  to  3  MeV  was  used. 
These  energies  cause  charge  penetration  to  centimeter  depths  in  the  test  materials  and  under  certain 
conditions  lead  to  tree  channel  breakdown  paths.  Advantages  of  this  method  include: 

i)  control  of  amount  and  spatial  distribution  of  charge 

ii)  direct  measurement  of  currents  during  propagatiopn 

iii)  direct  measurement  of  clannel  formation  velocities 

iv)  measurement  of  chaimel  morphology 

v)  production  of  extended  channels  for  comparison  to  models 

The  typical  procedure  for  implanting  charge  in  the  test  specimens  by  means  of  an  electron 
beam  has  been  described  in  detail  in  reference  [1].  These  smdies  employed  the  same  basic  methods 
with  a  beam  of  several  microamperes  ususally  directed  at  the  broad  face  of  a  sample  for  a  time 
sufficent  to  accumulate  a  net  charge  on  the  order  of  1  }iC/cm^,  Figure  1.  The  magnimde  and  spatial 
distribution  were  controlled  by  a  combination  of  beam  intensity  and  by  masking.  Discharges  could 
then  be  initiated  firom  various  locations,  for  example  from  one  end  as  in  Figure  2. 
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ELECTRON  BEAM 

lU'llllUUll  , 

1 _ 

7,89mm  thick 


Figure  1.  Sample  Geometry.  5cm  x  5cm  Width.  30mm  Dia.  Metallization 


charge  ao  iaplaatad  charge 


Figure  2  Electrical  "tree”  discharge  stops  propagation 
at  implanted  charge  boundary 
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2.5  New  Diagnostics 

In  order  to  better  determine  the  actual  distribudon  of  accumulated  charges  prior  to 
breakdown  it  was  necessary  to  employ  new  diagnostic  methods.  While  published  articles  had 
indicated  charge  accumulation  by  dosimetry  and  collected  current  methods  [2,3],  there  was  a  lack  of 
hard  evidence  for  true  charge  distributions  created  by  e-bcam  irradiation  without  insened  probe 
electrodes.  To  more  directly  measure  the  charges,  two  non-destructive  methods  were  used;  one  used 
the  Kerr  electro-optic  effect  and  the  other  a  newly  developed  electro-acoustic  technique.  Both  were 
successful  and  arc  described  further  in  Section  6. 

3.  CHARACTERISTICS  OF  SPACE-CHARGE-INDUCED  BREAKDOWN 

In  this  section  the  attributes  of  space  charge  induced  breakdo’vn  are  detailed.  The  branched 
breakdown  morphology,  its  depth  in  the  dielectric,  and  the  channel  sizes  and  lengths  are  presented. 
Then  events  during  channel  development  are  discussed,  such  as:  emitted  light,  currents  and  collected 
charge.  Finally,  material  specific  differences  are  discussed. 

3.1  Branched  Breakdown  Morphology 

The  dielectric  breakdown  induced  by  accumulated  space  charges  has,  in  every  case,  lead  to 
a  breakdown  path  which  was  an  extended  branched  structure.  A  permanent  damage  path  composed  of 
a  tree-like  set  of  interconnected  channels  is  produced.  Figure  2.  These  channels  do  not  recioss  the 
same  region  or  reconnect  one  another  to  form  loops.  An  exception  occurs  when  there  has  been 
multiple  charge  and  discharge  events  one  after  the  other  with  repeated  charge  implantation  into  the 
same  volume. 

3.1.1  Breakdown  channel  position 

Except  for  a  single  exit  path  to  the  surface,  the  physical  location  of  the  breakdown  channels 
resulting  from  e-beam  implanted  charges  was  normally  within  the  dielectric  and  mainly  within  the 
region  of  highest  space  charge  density.  The  location  of  the  channel  exit  point  and  the  breakdown 
initiation  was  often  stimulated  by  a  point  electrode  pushed  into  the  surface  of  a  charged  sample,  as 
illustrated  in  Figure  3.  Thfe  depth  below  the  dielectric  surface  where  the  channelling  occurred  was 
reasonably  confined  and  determined  by  the  energy  of  the  electron  beam.  Figure  4  shows  typical 
channel  depth  values  vs  e-beam  energy.  In  addition,  for  comparison,  expected  electron  range  values 
[1]  are  also  shown.  Note  that  the  channels  are  typically  at  a  depth  about  half  that  for  the  range.  This 
means  that  on  average  the  net  electronic  charge  accumulates  at  about  the  half  range  position  in  charge 
storing  dielectrics.  This  is  in  general  agreement  with  the  published  findings  of  Gross  and  Wright  [2] 
and  of  the  electro-accoustic  results.  Section  4. 
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Figure  3  Schematic  view  of  breakdovm  initiation 
and  propagation 
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3.12  Breakdown  channel  size  and  extent ' 

The  breakdown  channels  were  found  to  propagate  to  a  self-stopping  distance  which 
depended  upon  the  space  charge  density  in  the  material.  While  an  earlier  smdy  had  shown  the 
propagation  length  increased  with  electron  irradiation  dose  [1],  in  this  work  the  role  of  the  charge 
density  was  more  explicitly  identified.  In  particular,  the  charge  density  was  changed  in  one  set  of 
samples  by  means  of  cooling  the  samples  to  dry  ice  temperature.  -78*C.  The  cold  samples  had  greatly 
reduced  charge  mobilities  so  that  for  a  given  incident  charge,  the  cold  samples  contained  more 
accumulated  charge  and  as  a  result  exhibited  longer  channel  propagation.  This  difference  is  shown  in 
the  data  of  Figure  5. 

The  greater  charge  density  sustained  in  the  cold  samples  was  quantified  by  means  of  the 
electro-acoustic  charge  measurement  method,  as  described  in  Section  4  concerning  charge 
accumulation.  Samples  of  PMMA  were  measured  for  total  stored  charge  when  they  were  held  at 
various  temperatures.  The  charge  content  vs  time  for  these  cases  showed  that  at  room  temperature  a 
significant  amount  of  charge  was  already  transferred  out  of  the  sample  within  a  few  minutes  from 
e-beam  injection.  This  time  corresponds  to  the  amount  of  time  needed  to  perform  the  propagation 
length  test.  The  amount  of  charge  loss  is  about  25%  which  accounts  for  the  reduced  propagation 
length  at  room  teraperanire  shown  in  Figure  5. 

For  the  case  of  polycarbonate  material  the  results  were  differenn  In  a  manner  similar  to 
PMMA.  a  minimum  accumulated  charge  was  needed  to  initiate  breakdown  channels  in  polycarbonate. 
But  once  started,  they  propagated  for  very  long  distances,  even  1/2  meter,  at  moderate  charge 
densities.  The  difference  appears  to  be  the  result  of  a  significantly  lower  channel  impedance  caused  by 
carbonization.  Other  parameters,  such  as  currents  and  velocities  were  not  so  different  between  PMMA 
and  polycarbonate  as  will  be  discussed  later.  The  role  of  the  channel  impedance  is  discussed  in  the 
section  on  modelling,  Section  5. 

The  length  of  the  discharge  was  also  found  to  depend  slightly  on  the  width  of  the  charged 
region.  This  effect  is  shown  in  Figure  6.  Note  that  once  a  sample  was  4  to  5  cm  wide  there  was  little 
difference  from  that  for  a  much  wider  sample  of  10  cm  or  so.  The  reduction  for  narrow  widths 
appears  to  be  another  example  of  how  reduced  internal  fields  reduce  the  channel  propagation. 

3.1.3  Relation  to  fractals 

The  visible  branched  structure  of  the  breakdown  channels  looks  much  like  a  tree  or  an  aerial 
view  of  drainage  erosion  patterns.  These  structures  have  been  analyzed  in  mathematical  form  by 
means  of  a  formalism  called  fractals  as  developed  by  Mandelbrot.  [4]  For  the  case  of  space  charge 
breakdown  patterns,  three  regions  were  analyzed.  One  was  the  extreme  end  region 
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where  the  channels  stopped  propagation  even  thou^i  additional  space  charge  was  present.  Section 
3.1.2,  the  next  region  was  in  the  central  portion  and  the  third  near  the  discharge  root 

A  quantity  called  the  'branching  order'  is  used  to  determine  the  pattern  of  the  branching 
process.  For  the  case  of  the  tip  region,  an  expanded  view  of  the  channel  structure  is  shown  in  Figure 
7.  The  measured  branching  order  vs  position  is  shown  in  Figure  8.  The  near  constant  slope  indicates 
that  the  branching  process  is  a  regular  event. 

According  to  fractal  theory  the  negative  slope  corresponds  to  the  system  dimension.  One 
implication  is  that  similar  dimensions  correspond  to  systems  which  operate  under  the  same  basic 
physical  principal,  in  tl'is  case  that  energy  loss  is  minimized.  The  conclusion  is  that  the  fractal 
interpretation  lends  suppon  to  the  model  that  the  branched  breakdown  structure,  is  a  regular  pattern 
which  may  involve  an  energy  minimization. 

3.2  Processes  During  Channeling 

In  this  section  the  characteristics  of  events  which  occur  during  the  branched  channel 
propagation  are  presented.  Light  emission,  electric  currents,  charge  transfer,  and  the  role  of  local 
variations  in  charge  density  are  topics  that  are  presented. 

3.2. 1  Integrated  emitted  light 

To  help  establish  the  characteristics  of  the  channeling  processes  the  integrated  emitted  light 
produced  by  the  channeling  was  recorded  and  compared  to  the  physical  damage  caused  by  the 
channeling.  Specifically,  open  shutter  photographs  of  several  PMMA  samples  were  taken  during  the 
channel  development.  These  emitted  light  images  were  then  compared  to  subseqent  images  of  the 
physical  damage.  Figures  9  and  10  show  two  cases  of  emitted  light  and  physical  damage.  There  is 
very  good  correlation  between  the  two.  The  unusual  path  of  the  discharge  pattern  in  Figure  10  was 
made  by  masking  a  portion  of  the  sample  with  two  lead  disks  so  that  the  space  charge  was  excluded 
from  these  two  circular  regions. 

The  current  and  light  brighmess  appear  to  be  quantitatively  correlated  so  that  emitted  light 
may  be  used  as  a  good  diagnostic  for  discharge  damage.  Generally,  it  can  be  seen  from  the  figures 
that  the  brighter  regions  correspond  to  ths  'trunk'  charmels  where  larger  currents  of  longer  duration  are 
anticipated.  To  quantify  this,  induction  pickup  current  detectors  at  different  locations  along  a  sample 
should  be  compared  to  the  emitted  light  in  the  corresponding  regions.  Such  tests  were  beyond  the 
scope  of  this  program. 


2. 5 /lev 
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Figure  8  Branching  Analysis 


(a)  Open  shutter  emitted  light 
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(b)  Damage  Channeling 


Figure  9  Channeling  Process,  Uniform  Charge 
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3.2.2  Currents  during  channeling 

Tne  total  current  during  space  charge  induced  breakdown  events  was  recorded  and  showed 
a  remarkable  correlation  to  channel  development.  The  current  through  the  main  channel  and  out  the 
'exit'  path  was  measured  by  means  of  a  fast  pulse  Rogowski  coil  (10  ns  risetime)  around  a  shon 
length  conductor  attached  to  the  exit  path  (root)  region.  A  strip-line  geometry  was  used  to  keep  fast 
time,  nanosecond,  resolution.  Typical  currents  were  in  the  range  of  10  to  100  amperes  and  lasted  for 
about  200  to  800  ns.  For  a  uniform  implanted  charge  density,  the  current  immediately  rose  to  a  peak 
and  slowly  decayed  during  the  propagation.  Figure  11.  The  decay  of  current  corresponded  to 
inaeased  channel  length. 


3.2.3  Local  charge  density  variation 

More  importantly,  the  cunent  waveforms  were  very  consistent  with  the  local  charge 
density.  Thus,  when  the  implanted  charge  was  enhanced  or  reduced  in  a  local  region,  the  change  was 
observed  as  a  corresponding  change  in  the  current  waveform  at  the  channel  base,  ie.  root,  when  the 
discharge  dp  passed  through  the  region  of  different  charge  density.  The  synchronism  of  increase  and 
decrease  in  current  with  a  pattern  of  high  and  low  charge  density  regions  respectively  was  very 


strong.  Figures  12  and  13.  Because  of  the  rapid  response  in  root  current,  even  down  to  10  ns,  it 
appears  that  the  change  in  discharge  intensity  according  to  tip  region  conditions  is  evident  along  the 
whole  length  of  channel  connecting  to  that  region.  The  strong  control  of  discharge  propagation  by 
local  charge  variations  was  clearly  established  in  a  series  of  five  types  of  measurements; 

a)  step  increase,  Figure  12 

b)  4  step  ladder.  Figure  13 

c)  gap  jumpmg,  Figure  14 

d)  reservoir.  Figure  15 

e)  steering.  Figure  16. 

The  first  two  (a  &  b)  show  current  jumps  in  accordance  with  an  oscillating  in  space  distribution  of 
volume  space  charges.  The  next  (c)  shows  the  strong  reduction  in  current  as  the  propagation 
penetrates  a  short  region  of  no  charge.  The  sudden  jump  increase  in  current  when  a  discharge 
propagates  into  a  reservoir  of  charge  is  demonstrated  in  case  d).  And  finally,  the  path  of  propagation 
(e)  is  shown  to  accurately  follow  a  meander  path  of  implanted  charge. 

Overall,  there  is  very  consistent  evidence  that  the  local  charge  density  determines  the 
intensity,  direction,  and  (as  is  shown  in  the  next  section)  the  velocity  of  the  channel  propagation.  This 
finding  is  consistent  with  an  available  energy  model  for  channel  development  as  discussed  in  Section 
5. 
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Figure  11  Channel  propagated  full  length  of  lA  cm 
in  about  a  time  of  500  ns  according  to 
the  total  collected  current  waveform 
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Figure  16  Electrical  tree  follows  contour  of 
implanted  charge 


(r\\T-  HVRU 


Space  Charge  Report,  AFOSR 


Page  23 


3.2.4  Velocity  of  propagation 

The  velocity  of  propagation  of  the  channeling  discharge  was  measured  by  means  of  a 
relative  technique  which  compares  the  simultaneous  discharges  in  two  adjacent  paths,  the  so-called 
interferometer  method.  [1]  The  results  indicate  that  the  velocity  increases  as  about  the  square  of  the 
charge  density.  These  tests  were  made  with  cooled  samples  so  as  to  avoid  the  problem  of  charge 
leakage  during  setup.  Note  that  Polycarbonate  and  PMMA  exhibit  similar  results,  Figure  17. 

A  second  method  to  determine  velocity  is  by  means  of  the  'ladder'  pattern  of  Figures  13. 
Because  the  charge  spacing  is  at  equal  intervals,  the  time  between  current  peaks  is  inversely 
propordonal  to  velocity.  This  assumes  that  the  current  change  at  the  root,  where  the  current  monitor  is 
located,  does  not  have  a  delay  which  varies  appreciable  with  the  length  of  the  channel.  The  figures 
indicate  that  the  velocity  is  relatively  constant  with  channel  length  and  is  about  10^  cm/sec,  although  it 
docs  slow  down  somewhat  for  extended  lengths. 

3.2.5  Transferred  charges 

In  addition  to  tests  with  a  current  monitor  to  measure  the  total  current  to  the  root  of  the 
discharge,  some  tests  were  made  with  a  10  |ifd  capacitor  in  series  with  the  lead  so  as  to  measure  the 
net  charge  transferred  by  the  channeling  discharge  in  PMivLA.  The  results  of  these  tests  are  given  in 
Figure  18.  For  a  fixed  level  of  implanted  charge  density  the  amount  of  charge  transferred  is  relatively 
high  for  short  samples  and  becomes  progressively  lower  for  longer  samples.  Note  that  this  data  uses 
an  implanted  stored  charge  value  of  0.34  |iC/cm^  or  about  28%  of  the  incident  beam  for  the  irradiation 
level  of  3  kGys  used  on  these  samples. 

The  reduction  in  charge  transfer  efficiency  with  increased  propagation  length  is  expected 
from  the  fact  that  the  channeling  reaches  the  limit  of  its  extent  while  at  the  same  time  maintaining  a 
substantial  velocity  until  it  suddenly  stops.  That  is,  the  channeling  appears  to  be  less  effective  in 
gathering  energy  and  hence  charge  at  extended  distances.  The  reduced  charge  collection  is  also 
consistent  with  the  reduced  current  as  channels  become  long.  The  conclusion  is  that  a  substantial 
portion  of  net  change  remains  in  the  sample  in  regions  where  long  channels  propagate,  and  the 
efficiency  of  obtaining  electrostatic  energy  from  the  space  charge  is  reduced  by  long  channels  in 
PMMA. 


The  charge  collection  efficiency  for  fixed  size  disk  samples  of  46  cm^  irradiated  active  area 
was  measured  vs  the  total  electron  injection.  The  results  are  shown  in  Figure  19.  Note  that  some  of 
the  reduction  in  the  fraction  of  charge  collected  at  higher  injected  charge  levels  is  the  result  of  less 
fractional  stored  charge  at  the  moment  of  measurement,  see  sections  3.2.6  and  4.1.2. 
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Figure  17  Channel  Development  Propagation  Velocity 


Xrracilay  (om') 

I 

Figure  18  Transferred  Charge  in  PMM.A.  j 

(^/t-WKRL  i 


j 

1 


'_N  .*  j< 


x”. 


Page  >.0 


Z5/^eV 

P/'IAlft;  ‘trhlc.k 

4 


0  J  O.H- 

0.8  l.l  J.Q> 

^o  zi 

HaiiMr  hm'it 

XwTEcTED  Charge 

Figure  19 

Charge  Collection  Efficiency  for  46  cm^ 
round  Charged  Area,  PMMA 

P^ir-HVRL 

_V7?_#:  /  ji...  r.....-V..  *.  fkJ?, . 

Space  Charge  Repon,  AFOSR 


Page  27 


3.2.6  Incident  flux  rate  effects 

By  modification  of  the  intensity  of  the  e-beam  it  was  possible  to  change  the  rate  at  which 
charge  was  implanted  in  the  experimental  samples.  The  electro-acoustic  method  was  used  to  examine 
the  effect  on  the  internal  charge  accumulation.  Figure  20  shows  the  measured  internal  charge  vs  the 
total  incident  e-bcam  charge.  Four  different  intensities  were  used  ( as  indicated  by  the  numbers  beside 
the  data  points.  At  low  incident  charge  levels,  much  of  the  charge  is  crapped  and  there  is  a  linear 
increase  of  stored  charge  with  trapped  internal  charge.  However,  at  higher  incident  charges  less  of  the 
charge  is  accumulated. 

The  lower  efficiency  of  charge  trapping  at  the  higher  implant  intensities  can  be  expected 
from  the  characteristics  of  charge  transpon.  At  the  higher  intensities  a  greater  amount  of  power  is 
delivered  to  the  sample  and  heating  occurs.  But  there  is  a  strong  dependence  of  charge  transpon  on 
temperature,  as  shown  in  Section  4.1.2,  so  that  during  the  charge  implant  time  a  significant  amount  of 
charge  is  conducted  from  the  volume  by  transpon.  A  second  means  by  which  reduced  net  charge  is 
accumuniated  at  high  implant  rates  is  due  to  radiation  induced  conductivity.  The  usual  model  for  this 
effect  is  that  this  form  of  conductivity  is  proportional  to  radiation  dose  rate.  Hence,  at  high  implant 
rates  a  greater  percentage  of  the  injected  charge  is  conducted  from  the  volume. 

3.3  Material  Effects 

In  this  section  the  influence  of  the  dielectric  material  on  the  channel  development  is 

discussed. 

3.3.1  Initiation,  channel  propagation 

There  is  a  substantial  difference  in  the  amount  and  morphology  of  the  branched  discharge 
channels  which  develop  in  space  charged  regions  according  to  the  dielectric  materiaL  These  differences 
will  be  discussed,  but  first  it  is  worth  pointing  out  that  the  basic  structure  of  this  class  of  discharges  is 
remarkably  the  same,  independent  of  material.  Specifically  all  materials  studier*  show  common 
fundamental  characteristics  of: 

-  the  sudden  formation  of  channel  discharges  at  sufficient  charge  levels 

-  greater  ability  to  accumulate  charge  at  lower  temperatures 

-  light  emiaed  during  channel  development 

-  similar  velocities  of  propagation 

-  discharge  channels  which  do  not  recross  a  region  during  one  event 

-  triggering  of  discharge  by  local  high  elecaic  fields 

-  e-beam  caused  space  charge  acctunularion 

-  increased  depth  of  channels  below  surface  with  greater  implant  energy 
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Figure flU)  Initial  Trapped  Charge  from  1  MeV  Electrons  in  PfiMA 


Page  29 

The  differences’  associated  with  materials  are  most  notable  in  three  categories: 

•  iaidadon  condidon 

-  propagadon  extent 

-  channel  damage  morphology 

The  cause  for  these  differences  seems  to  be  associated  with  two  processes:  the  charge  accumulation 
physics  and  the  channel  plasma  physics  and  damage.  The  charge  accumulation  at  room  temperature  in 
materials  such  as  epoxy  can  be  very  small  compared  to  PMMA.  However,  with  epoxy  it  was  possible 
to  obtain  sufficient  stored  charges  to  cause  channel  discharges  by  first  cooling  the  sample  to  -78’C. 
The  difference  in  charge  storage  accounts  for  the  initiation  and  some  of  the  propagation  differences 
among  the  different  materials. 

The  propagation  extent  is  determined  not  orJy  by  the  amount  of  space  charge  available,  but 
also  apparently  by  the  conductivity  of  the  channel  during  propagation.  For  example,  polycarbonate 
(LEXAN)  forms  black  (carbon  looking)  charmels  and  this  material  exhibits  unsually  long  channels, 
even  for  charge  levels  slightly  above  threshold  values.  This  is  interpreted  to  mean  that  the  channel 
impedance  is  low  even  for  long  channels  and  hence  discharge  propagation  can  continue  because  the 
available  energy  remains  high.  Some  channels  propagated  the  full  length  of  a  76  cm  long  sample 
with  only  0.4  4C/cm^  implanted  charge.  For  comparison,  the  same  implanted  charge  in  PMMA  yields 
a  channel  propagation  of  just  3.6  cm  when  it  self-stopped. 

Surprizingly,  even  though  the  propagation  length  is  very  different,  the  velocity  of 
propagation  in  polycarbonate  is  similar  to  that  in  PMMA.  Tests  with  periodically  spaced  positions  of 
high  charge  density,  as  in  Section  32.3,  showed  that  the  velocities  and  currents  in  these  materials 
were  different.  Figure  21,  with  the  larger  values  in  LEXAN.  The  higher  currents  confirm  the  earlier 
statement  of  lower  channel  impedance  in  this  material.  The  relative  velocity  measurements  with  the 
'interferometer'  struemre.  Section  3.2.4,  showed  a  very  similar  dependence  of  velocity  on  implanted 
charge  level,  Figure  16. 

Not  all  insulators  readily  accumulate  charges.  Polyethylene,  for  example,  is  a  relatively 
good  dielectric.  Yet  at  room  temperature  it  did  n2l  form  discharge  channels  for  e-beam  intensity  values 
similar  to  those  used  in  the  PMMA  or  polycarbonate  tests.  Electro-acoustic  measurements  showed  that 
the  stored  charge  in  polyethylene  was  indeed  orders  of  magnimde  below  that  of  the  other  two 
materials.  Also  the  charge  decay  rate  was  much  higher.  Tests  were  also  made  with  Mylar  polyester. 
However,  it  was  only  available  in  thin  sheets  so  that  even  at  200  keV  beam  energy,  most  of  the 
electrons  penetrated  the  sample.  Hence,  there  was  limited  charge  accumulation  and  internal  discharges 
did  not  occur. 
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Figure  21  Channel  Propagation  vs  Material  6  gap 

oscillating  charge  of  0.6  yC/cm^,  2.5  MeV 
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3.4  Embedments 

3.4.1  Spheres 

A  sample  was  with  small  3.2  mm  diameter  metal  spheres  embedded  in  a  plane  at  about 
the  expected  depth  of  the  subsequendy  implanted  charge,  ie.  about  8  to  9  mm  below  the  injection 
surface.  Figure  22  shows  the  channel  structure  and  the  termination  of  each  channel  which  strikes  a 
sphere.  This  channeling  termination  is  taken  as  evidence  for  a  greatly  reduced  electric  tip  field  due  to 
the  larger  radius  of  curvature  of  the  sphere  compared  to  the  natural  channel  tips. 


4.  MODEL  FOR  CHARGING  OF  DIELECTRICS  BY  E-BEAM  IRRADIATION 

In  this  Section  the  details  of  a  model  for  the  charge  accumulation  process  in  solid  dielectrics 
resulting  from  controlled  charge  implantation  by  e-beam  irradiation  are  presented  along  with 
experimental  results  which  support  the  models  developed. 

4.1  Internal  Charge  Decay 

Assuming  that  a  net  charge  accumulation  within  a  solid  dielectric,  such  as  caused  by  e-beam 
irradiation,  has  occurred,  the  charg  will  subsequently  move  and  transpon  out  from  the  volume.  In 
dus  section  the  net  charge  decay  within  the  dielectric  is  modelled.  This  is  analysis  is  made  in  order  to 
understand  the  net  charge  which  develops  during  the  charging  process  as  well  as  to  understand  the 
subsequent  time  dependence  of  changes  in  volume  charges. 

4.1.1  Decay  with  time,  initial  charge  effects 

For  various  initial  charge  levels  within  the  dielectric,  as  determined  by  the  e-beam 
implantation,  the  decay  of  net  volume  charge  in  samples  held  at  room  temperature  was  measured  over 
a  period  extending  to  two  months.  Figure  23.  The  samples  remained  shorted  on  all  sides  by  wrapping 
with  metal  foU  during  the  storage  interval  between  measurements.  Two  characteristics  are  readily  seen 
in  these  curves.  First,  the  relaxation  time  grcalty  depends  upon  incident  beam  charge.  For  example, 
whereas  the  sample  irradiated  to  (  3  kGy )  decayed  its  charge  by  a  factor  of  50  in  1,000  hours,  the 
sample  irradiated  with  about  (  70  Gy  )  showed  less  than  a  factor  of  2  charge  decay  in  the  same 
1,000  hours.  Second,  all  curves  appear  to  extrapolate  to  a  single  asymptote  of  the  form  1/  St,  where 
sis  a  constant  1.6  x  10^  (C-hr)'^  and  t  is  time  in  hours  for  the  7  cm~  area. 

4.1.2  Temperature  effects  an  decay 

Additional  samples  stored  at  60*C  and  at  90’C  were  measured  for  charge  relaxation  with  the 
same  implanted  charge  of  (3  kGy  ),  as  that  for  the  reference  room  temperature  test  (~  28*  C).  In 
addition,  a  sample  at  dry  ice  temperature,  -80*C,  was  irradiated  to  ( 1.5kGy ).  The  resultant  decay 
curves  are  shown  in  Figure  24.  Solid  lines  are  drawn  with  a  shift  of  a  factor  of  10  for  each 
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Figure  22  Discharge  Channels  in  PMMA  Semple  with  3.2mm  Die.  Metal  Sphere 

Embedments.  Showing  Channels  Terminate  at  Spheres 
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Figure  23  Charge  Relaxation  in  PMMA  after  Irradiation  at  1  MeV 
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Figure  U.  Acoustic  Signal  Waveforns  for  Cliarge  I.aycr  In  I'lOIA 
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SPACE  CHARGE  ACCUMULATION  IN  POLTMETHILMETHACRYLATE 
RESULTING  FROM  ELECTRON  BEAM  IRRADIATION 

C.M.  CooROk  M.  Zahtt*  E.A.  Wright,  and  M.  Hihita 
Maasachusatts  Inatitnte  of  Technology 
Laboratory  for  Electronagnetio  and  Electronic  Systems 
High  Voltage  Research  Laboratory 
Cambridge.  MA  02139 

INTRODUCTION 

The  measnrement  of  volume  charges  within  insulators  provides  needed 
information  important  to  spacecraft  application.  The  actual  charge  values 
help  in  establishing  models  for  charge  accumulation  and  decay  and  lead  to  a 
better  understanding  of  spontaneous  discharges  of  insulators  exposed  to 
energetic  electron  irradiation.  In  this  work,  two  methods  of  measnrement  of 
trapped  charges  in  polymethylmethacrylate  are  used;  the  electro'-aeoustie 
method  [1]  and  the  Kerr  electro-optic  method  [2].  Both  provide  quantitative 
information  about  internal  space  charges. 

Of  particular  interest  are  the  volume  charges  in  PMMA  resulting  from 
exposure  to  megavolt  energy  electrons  [3.4.5.61.  A  series  of  tests  were  made 
to  determine  how  charges  accumulate  according  to  dose,  what  are  typical 
profiles  of  charge  with  depth  into  the  bulk  material,  and  what  are  typical 
relaxation  decay  times  for  this  charge  to  leave  the  sample. 

EXPERIMENTAL  METHOD 

The  electro-acoustic  method  yields  volume  charge  information  by  applying 
a  short  pulse  electric  field  to  the  sample  and  measuring  the  resultant 
mechanical  pressure  wave  which  arrives  at  the  electrodes  [1.7].  The  shock 
magnitude  and  time  delay  of  arrival  are  measured  with  a  piezoelectric 
transducer  and  provide  information  about  the  size  and  distribution  of  the 
volume  charge.  Figure  1  shows  a  schematic  of  this  type  of  system. 


The  Kerr  method  employs  incident  light  polarized  at  45°  to  the  direction 
of  the  electric  field  in  the  sample  and  an  analyzing  polarizer  placed  after 
the  sample.  [2,5]  The  resnltant  transmitted  light  intensity  is  a  function  of 
the  internal  field  and  when  uniformly  illuminated  provides  a  complete  field 
intensity  map  over  the  depth  of  the  saaiple.  Figure  2  shows  a  schematic  of 
apparatus  for  the  Kerr  measurements. 

Both  methods  as  a  practical<matter  need  a  calibration  constant.  This  is 
determined  by  means  of  tests  with  a  known  applied  field  [2,7].  With  the 
separately  determined  constants,  the  measurements  become  quantitative. 

The  measurements  were  made  on  PMMA  samples  typically  0.6  to  1.3cm  thick. 
Electrons  from  a  Van  de  Graaff  generator  entered  the  saaqtle  on  the  broad  area 
face  with  energies  in  the  0.5  to  2.5MeV  range  with  current  densities  5  to 
250nA/cm*.  The  samples  had  thin  metal  foil  electrodes  which  short-circuited 
top  and  bottom  faces  throughout  the  tests. 

VOLUME  CHARGE  ACCOMUUTION 

The  net  accumulation  of  charge  within  the  samples  when  at  room 
temperature  was  determined  by  the  electro-acoustic  method  and  was  found  to 
increase  linearly  with  dose  up  to  about  0.1  or  0.2  kGy.  At  higher  doses,  a 
saturation-type  response  occurs.  Figure  3.  This  may  result  from  one  or  more 
of  the  following  processes:  space  charge  repulsion  by  the  already  trapped 
volume  charge,  radiation  induced  conductivity,  snd/cr  heating  which  reduces 
relaxation  times, 

VOLUME  CHARGE  PROFILE 

The  electro-optical  method  was  used  to  reveal  the  distribution  of  volume 
charges  by  virtue  of  the  electric  fields  which  they  induce.  Figure  4  shows  a 
typical  measured  field  distribution  for  2NeV  electrons  incident  on  a  ll.lwa 
thick  sample.  The  influence  of  incident  beam  energy  is  shown  in  Figure  5, 


where  the  depth  for  zero  field,  i.e.,  the  point  of  internal  voltage  maziana, 
is  plotted  against  beaa  energy.  A  linear  increase  of  depth  with  energy  occurs 
in  this  case  where  the  saaple  is  thicker  than  the  electron  range.  Note  the 
aeasnred  distribution  is  very  close  to  that  theoretically  predicted  by  Tanaka, 
et  al.  [8]  and  to  measureaents  obtained  with  collecting  plates  [9].  which  is 
good  confiraation  of  the  results. 

VOLUME  CHARGE  DECAY 

The  electro-acoustic  method  was  used  to  measure  the  remaining  net  charge 
as  a  function  of  time  after  a  given  short-tera  injection  of  charge.  Figure  6 
shows  the  results  for  the  case  of  decay  after  incidence  of  a  l.ONeV  beaa  to  a 
level  of  10~‘coul/ca* .  The  decay  to  half-value  occurred  in  about  25  to  30 
minutes,  but  decay  to  one-tenth  initial  value  took  about  70  hours.  Seduction 
by  a  factor  of  100  is  extrapolated  to  about  3.000  hours.  This  is  a  non¬ 
exponential  decay  and  at  longer  times  follows  the  self-precipitation  model  of 
a  (t+t)“^  time  dependence. 

CONCLUSIONS 

Volume  space  charge  accumulation  in  the  insulating  dielectric  PHMA  has 
been  shown  to  occur  as  a  result  of  electron  irradiation.  At  low  doses  the  net 
charge  increases  linearly  with  dose,  but  charge  decay  can  be  significant  and 
acts  to  continuously  discharge  the  saaple.  At  high  doses,  a  saturation  of  net 
charge  occurs  quickly.  At  low  and  moderate  doses,  the  charge  profile  in  depth 
appears  to  be  consistent  with  model  predictions. 
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Figure  5,  Depth  into  PMMA  for  Zero 
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Figure  5,  Volume  Charge  Decay  in  PMFIA  after 
Irradiation  with  1  MeV  Electrons  to 
3  k6y  in  lOsEc,  Electro-acoustic  i-Iethod 
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307C  change  in  temperature.  Note  the  1/  3t  asymptotes  are  shifted  by  a  factor  of  10,  too.  The  data  at 
dry  ice  temperature  showed  no  significant  decay  even  though  the  total  charge  of  2  |iC  ( 0.29  p.C/cm“) 
was  relatively  high. 


4.1.3  Mobility  model 

The  dependence  of  charge  relaxation  time  on  net  volume  charge  is  a  phenomena  which 
contradicts  traditional  constant  conductivity  models  used  to  describe  "resistive"  decay  with  a 
magnitude  independent  time  constant  of  ~  h .  The  charge  dependent  conductivity  data  suppons  a 
mobilitv  model  of  the  form: 


J  =  piiE 


where  2  is  the  volume  charge  density,  |j.  is  the  carrier  mobility,  and  E  is  the  electric  field.  Continuity 
and  Poissons  equation  for  a  single  carrier  species  yields: 


|£-  +  i  p>  +  .  Vp  =  0 


For  a  small  gradient  in  charge  density,  the  solution  is: 


P  “  — —  =  Po  TIT 


e 


t+T 


where  p  is  the  initial  charge  density  and  x  is  the  relaxation  time: 


II  A 

fro 


which  depends  upon  p .  For  the  simple  case  of  a  uniform  u  and  a  single  dominant  carrier 
integration  over  the  thickness  yields: 


,t+l 


(1) 
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where  q^s  the  inital  charge  per  unit  area  and  ®  is  a  constant: 

/a  ^  1 

a  •  —  (coul-sec)'^ 

'  sD 

The  equation  for  charge  relaxation,  Eq.  ( 1 ),  is  plotted  in  Figure  25  and  appears  to  be  consistent  with 
the  low  implanted  charge  data  of  Figure  20.  At  higher  implanted  charge  levels,  accelerated  charge 
relaxation,  perhaps  associated  with  radiation  induced  conductivity  or  multiple  carrier  species  appears  to 
dominate  initially,  and  after  about  1,000  hours  the  self  relaxation  at  constant  mobility  is  seen. 

The  thermal  effects  may  be  modelled  as  a  thermally  dependent  mobility,  where: 

with  T  in  degrees  absolute  and  A  and  Vg  constants.  For  a  factor  of  10  change  in  mobility  with  a  30*C 
change  near  room  temperature,  Vg  has  a  value  of  about  0.88  eV.  Hence,  from  the  room  temperature 
data,  p.  =  8.4  X  10‘^°  m^/v-sec.  .And  with  the  temperature  effect  included: 

4  =  8.2  X  10-9  jS  e  -Vg/TcT 


Funher  investigation  of  this  model  coula  be  made  especially  with  consideration  of  varying  the  rate  of 
implanting  charge.  Also  adjustments  for  the  case  of  high  total  injected  charge  could  be  smdied. 


To  show  the  efficiency  of  the  charge  crapping  process,  the  net  accumulation  of  charge 
within  the  samples  at  almost  constant  irradiation  time  =  30  sec,  when  at  room  temperature,  was 
measured  and  found  to  increase  linearly  with  beam  charge  (injected  charge)  up  to  about  0.1  kGy.  At 
higher  beam  charges,  a  saturation-type  response  occurs,  see  Figure  20.  This  may  result  from  one  or 
more  of  the  following  processes:  charge  repulsion  by  the  already  trapped  volume  charge,  radianon 
induced  conductivity,  and/or  heating  which  reduces  relaxation  time. 


Laver 


Typical  mean  charge  depths  were  measured  for  several  e-beam  implantation  energies.  The 
results  show  that  the  mean  charge  position  is  at  a  depth  somewhat  greater  than  the  channel  position  but 
less  than  the  electron  range  in  the  material,  see  Figure  4.  This  value  is  consistent  with  e-beam 
stopping  calculations  published  in  die  literature.[10] 
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5.  MODEL  FOR  BREAKDOWN  THROUGH  VOLUME  SPACE-CHARGED  REGIONS 


5.1  ErgpagariQn  Model 

5.1.1  Absolute  available  energy 

Because  there  is  no  applied  electric  field  from  an  external  power  supply,  the  discharges 
investigated  in  these  studies  demonstrate  very  clearly  that  energy  to  propagate  a  breakdown  event  can 
come  from  static  charges  distnbuted  in  a  volamc.  How  the  discharge  propagates  can  be  interpreted  as 
based  upon  electric  fields  and  the  forces  they  exen  on  charge  carriers  or  as  electrostatic  energy  which 
is  convened  to  light  and  heat  as  the  system  lowers  its  overall  internal  energy. 


In  the  field  driven  model,  a  discharge  propagates  when  the  field  resulting  from  the  trapped 
charges  exceeds  a  threshold  field  value,  E^j^,  necessary  to  compensate  for  loss  processes.  For  the 

self-breakdown  case,  space-charges  and  their  images  at  surfaces  are  the  sole  source  of  the  field.  This 
means,  for  a  uniform  initial  charge,  that  the  internal  field  depends  linearly  on  charge  density.  For  a 
uniform  layer  of  thickness  w,  remaining  length  H,  charge  densin  :  ,  as  depicted  in  Figure  26,  The 

space  charge  field.  Ej,  at  the  edge  is: 


"  ne 


thus : 


,2H.  z  ^  .2z. 

arcctg  -  —  arcctg  (— ) 


1,.  [,i  .  [(i.  (il'l) 


♦u; 


XIA  AilAO  LLl' 


is  model  the  channeling  should  continue  as  long  as  charge  was  available.  To  account  for  the 


experimental  finding  of  self-stopping  of  channels  in  materials  such  as  PMMA,  Section  3.1.2,  it  is 
necessary  to  have  a  means  oy  which  less  field  is  available  at  tr.e  propaganng  front  of  the  discharge. 
Two  mechanisms  for  this  ha\c  been  identined,  and  they  arc  discussed  further  in  the  next  section 
(5.1.2). 


For  the  energy  model,  it  is  assumed  that  the  power  to  sustain  the  discharge  and  its  losses 
[5,6]  is  derived  from  the  change  in  electrostatic  energy  before  and  after  the  breakdown  channeling  has 
past  through  a  region  times  the  velocity  of  propagation. 

=  AUjfz)  X  V£(z) 


Figure  26  Geometry  for  Edge  Field  Calculation 


,2cr>nzaArwrd 


Space  Charge  Report,  AFOSR 


Page  ^0 


where  P^.  is  the  power  needed  per  unit  propagation  in  the  z-dircction,  I  1^\^/  aS  uic  viwvuusiEtic  energy 
per  unit  distance  in  the  z-direcrion,  and  Vf  is  the  propagation  front  velocity. 


Hence  a  principle  question  is  the  amount  of  energy  available.  To  determine  this,  the 
electrostatic  energy  associated  with  an  implanted  charge  and  the  energy  after  a  channel  front  has 
propagated  through  the  charge  is  determined  so  as  to  evaluate  the  change  in  energy.  To  obtain  a  value 
for  the  energy  after  a  discharge  has  occurred,  a  model  for  the  channeling  was  used.  It  assumes  a 
conducting  plane  at  potential, taken  to  be  zero;  ie.  the  channel  impedance  is  zero.  A  conducting 

plane  is  used  as  an  equivalent  for  the  discharge  which  is  actually  composed  of  many  closely  space 
channels. 


A  calculation  for  tiie  change  in  electrostatic  energy  by  the  discharge  was  made  assuming  a 
geometry  similar  to  the  experiments.  A  dielectric  sheet  of  material  of  thickness  2Rq,  is  used  where  the 

distance,  Rq,  is  the  maximum  penecranon  distance  of  the  charge  layer.  The  charge  is  acmally  located 
mainly  in  a  layer  of  nnite  thickness,  but  for  the  calculation  the  charge  is  assumed  to  be  distributed  in 
one  of  two  forms:  constant  or  linearly  increasing  with  depth,  both  to  a  maximum  depth  of  Rq.  The 

configurations  and  results  are  shown  in  Figures  27  &  28  .  For  both  it  is  seen  that  the  energy  change 
for  discharge  channels  at  a  depth  a  from  the  surface,  is  sharply  at  a  maximum  for  specific  positions  of 
the  discharge,  x=a  /Rq.  Specifically,  it  is  seen  that  energetically  the  discharge  will  strongly  prefer  to 

be  at  a  depth  of  about  0.7  times  Rq.  This  value  is  in  close  agreement  with  the  experimental  findings. 
The  equations  for  the  initial  energy,  Uq,  and  the  change  in  energy,  AU,  are  for  the  case  of  uniform 
charge  density  to  Rq: 


Pn  -  (1  +  .5.) 


^  1-  fx  = 

24e  2  L 


+  + 


While  for  the  case  of  linearly  increasing  charge  density,  these  values  are: 


Figure  28  Calculated  Energ/  Preferred  Discnarge Channel  Location, 
Linearly  Increasing  Charge  Density  to  a  Depth  Equal  to 
the  Range  in  a  Material  of  Thickness  2  times  the  Range 
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5.1.2  Relative  available  energy,  Channel  impedance  effects 

While  the  above  model  describes  the  depth  of  the  channeling,  it  is  necessary  to  account  for 
the  self-stopping  length.  If  the  efficiency  of  charge  neutralizdon  by  the  discharge  channels  is  not 
unity,  then  a  remanant  charge  remains  and  there  is  a  lower  field  at  the  propagating  discharge  front. 
But  another  way  to  model  self-stopping  effects  is  to  allow  a  finite  impedance  to  the  discharge  channels 
so  that  there  is  a  voltage  drop  along  the  length  of  the  channels. 


Using  a  one-dimensional  model,  the  tip  potential,^  will  differ  from  the  ground  potential, 
0,  at  the  beginning  root  of  the  discharge.  It  will  increase  with  channel  length  according  to  the  integral 
of  voltage  drop  per  unit  length  of  channel,  Eg,  over  the  length  of  the  channel. 

0^=  Eg(l,t)  dt 


Expected  values  for  Eg,  according  to  the  literature,  [5,6]  indicate  that: 

where  j  is  the  channel  current  density.  Expressed  in  another  way  the  channel  resistance  per  unit  length 
varies  inversely  with  the  square  root  of  current  density.  Now  letting  the  channel  propagation  velocity 
be  proportional  to  the  available  energy  per  unit  length  of  travel,  we  have  that  velocity  is  initially 

proportional  to  the  charge  density  squared.  Tnis  is  because 

> 

AU  =  K  [  pj-  ] 

As  the  channel  becomes  long,  the  finite  resistance  lowers  the  tip  potential,  lessening  the  energy  per 
unit  path  avaHabie,  this  lowers  the  channel  current  density  at  the  channel  tip.  But,  because  it  was 
shown  by  the  variable  implanted  charge  density  studies  that  the  whole  length  of  channel  quickly 
responds  to  changes  at  the  tip,  the  whole  channel  will  also  have  a  reduced  current  density. 

Now  this  provides  rapid  feedback  to  cause  even  greater  tip  potentials  and  hence  lower 
energy  available  because  the  lower  current  density  raises  the  channel  resistance  and  thereby  raises  the 
voltage  drop.  This  progessivc  reduction  in  channeling  rapidly  comes  to  a  complete  halt  when  the 
available  energy  per  unit  length  of  growth  drops  below  the  threshold  power  needed  to  sustain  the 
channel  development.  Another  factor  which  may  occur  is  a  change  in  the  needed  power  to  sustain  a 
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a  channel  due  to  its  ablation  of  the  material  which  expands  the  channel  diameter. 

Thus,  channel  propagation  is  modelled  as  being  driven  by  the  available  energy  associated 
with  the  channel  developing  into  a  charged  region  exceeding  the  power  expended  for  that 
development.  Furthermore,  to  explain  the  experimental  findings  of  well  defined  self-stopping  lengths 
for  some  materials,  a  finite,  current  dependent  resistance  of  the  channels  is  employed. 


A  fundamental  long-term  failure  process  in  solid  dielctrics  often  occurs  as  an  internal 
branched  channel  structure,  which  because  of  its  shape  is  referred  to  as  an  'electrical  tree'.  This  tree 
channeling  usually  develops  over  a  period  of  time  at  stress  with  the  channel  branches  progressively 
extending  until  one  or  more  in  series  bridge  the  distance  bewteen  the  electrodes.  The  failure  is 
observed  as  electrical  conduction  via  a  connecting  path  through  one  or  more  of  the  interconnecting 
channels.  Cables  bushings,  capeitors  and  other  high  voltage  apparatus  can  fail  because  of  the  treeing 
process. 

The  extension  of  existing  branches  and  the  development  of  new  branches  in  an  electrical  tree 
results  in  the  tree  'growth'  or  expansion.  This  growth  is  proposed,  based  upon  the  results  of  these 
investigations,  to  occur  as  a  three-step  process  which  is  composed  of: 


1)  charge  injection  and  accumulation 

2)  dscharge  initiation 

3)  discharge  propagation 


The  first,  injection  and  accumulation,  is  dependent  upon  insulator  material,  electrodes,  applied  electric 
stress  and  temperanu'e.  The  second,  initiation,  is  dependent  upon  the  magnitude  of  the  first  stage  and 
on  the  applied  stress,  while  the  third,  propagation,  is  relatively  independent  of  temperature  and  the 
type  of  initiation.  The  three  steps  normally  repeat  rapidly  so  as  to  yield  an  almost  continuous  looking 
tree  growth  over  a  period  of  time.  Alternately,  the  growth  may  exhibit  distinguishable  steps  with  an 
apparent  pause  between  steps  associated  with  the  time  for  charge  accumulation  and  initiation. 


As  an  exiunple  of  electrical  treeing  and  a  comparison  with  channel  discharges  created  by 
e-beam  implanted  charge.  Figure  29  shows  side  views  for  trees  originating  from  high  field  point 
electrodes.  A  pointed  electrode  brought  adjacent  a  small  beam  charged  region  is  used  for  comparison. 


U) 


(b) 


Figure  29  Electrical  Discharge  Channels  in  PMMA,  (a)  AC  Voltage 
with  3mm  Point- to-Plane Gap,  (b)  s-beam  Implanted 
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6.  NEW  DIAGNOSTICS 

Two  methods  of  measurement  of  implanted  charges  in  solid  dielectrics  were  used  to 
evaluate  the  charge  magnitude  and  distribution.  One  employed  the  Kerr  electro-optic  effect ,  ie. 
birefringence,  while  the  other  employed  a  new  ciectro-acoustic  method.  Both  arc  described  below. 

6.1  Electro-optic  Measurements  of  e-.beam.  Implanted  Charges 


The  electro-optic  method  yields  volume  charge  information  by  applying  incident  light 
polarized  preferably  45*  to  the  direction  of  the  internal  electric  fields  in  the  sample.  After  passing 
through  the  sample,  the  light  is  passed  through  a  second  polarizer  at  90*  to  the  first.  This 
configuration  is  depicted  in  Figure  30.  [7]  'Fhe  resultant  trairsmittcd  light  intensity  is  modulated  by  the 
internal  field  so  that  by  examining  the  fringe  pattern  a  complete  mapping  of  the  internal  fields  is 
obtained.  Some  difficulty  was  encountered  with  the  mechanical  birefringence  which  introduced  an 
initial  fringe  pattern.  Simple  thermal  annealing  did  not  remove  these  mechanical  stresses  so  that  it  was 
necessary  to  subtract  their  effects  by  assuming  that  they  remained  constant.  Monchromatic  green  light 
from  a  mercury  discharge  lamp  was  used  to  remove  the  light  wavelength  dependence  of  the 
birefringence. 


Typical  results  are  shown  in  the  photographs  of  Figure  31  where  a  portion  of  the  sample 
surface  was  irradiateH,  Tlie  changing  patterns  during  the  decay  period  after  irradiation  is  stopped, 
clearly  show  the  initial  field  and  its  relaxation  as  the  charges  leave,  Figure  32,  in  addition  to  the  initial 
mechanical  stress.  The  dark  and  light  regions  of  the  fiinge  pattern,  min  and  max,  show  the  field 
exceeds  a  MV/cm.[7]  The  line  of  zero  field,  ie.  the  potential  maximum,  lies  below  the  surface  as  does 
the  charge  max  at  the  fiinge  density  maximum,  in  correspondence  to  the  discharge  channeling  depth. 

6.2  Electro-acoustic  Method 

The  electro-acoustic  method  was  developed  by  Prof.  Takada  from  Musashi  Institute  of 
Technolof^'  in  japan.  [  8]  He  came  to  MIT  for  a  technical  visit  and  collaborated  in  making  test 
measuiments  for  this  study  on  samples  implanted  with  electrons  by  the  MTT  accelerator. 


The  method  yields  volume  charge  information  by  applying  a  short  pulse,  typically  about  90 
ns,  electric  field  to  the  sample  and  measuring  the  resultant  mechanical  pressure  wave  which  arrives  at 
the  electrodes  [8,9].  The  shock  magnitude  and  time  d<5iay  of  arrival  are  measured  with  a  piezoelectric 
tranducer  and  provide  information  about  the  size  and  distribution  of  the  volume  charge.  Figure  33 
shows  a  schematic  of  the  system.  A  typical  acoustic  pressure  waveform  for  an  e-bcam  irradiated 
sample  shows  the  three  pulses  teo  for  the  surface  image  charges  and  one  for  the  implanted  charge  peak 
in  the  volume,  Figure  34.  Two  publications  about  this  method  have  been  made  and  are  attached. 


Alii  or 


Fijuro  13J  A  jtrejsed  photoelastic  model  in  a  plane  polariscope. 
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Ff|iir*  13.7  A  stressed  photoelastic  model  in  a  drcular 
and  analyzer,  crossed  quarter-wave  plates). 
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Figure  30  Electro-Op  lie  Method  of  internal  Charge  Measurement 
(from  Optical  Methods) 
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Figure  3 1  Localized  Irradiation  of  PliMA  Sample,  1 2.7mm  Thick 
and  1 0cm  Long  Optical  Path  Showing  Fringe  Shift 


tr\(L 


Figure  32  E-beem  Implanted  Charge  Deco/  in  PMMA,  Visualized  by  Optical  Method 
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Figure  34  Waveform  from  Pulse  Electro-Acoustic  Char  ge  Detection  System 
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7.  CONCLUSIONS 


The  investigations  were  successful  and  yielded  new  information  about  fundamental  physical 
processes  responsible  for  major  limitations  in  the  performance  of  electrical  insulation,  especially  for 
high  stress  applications  such  as  used  in  pulse  power  and  space  satellites.  New  models  to  acccount  for 
charge  accumulation  in  the  volume  of  dielectrics  were  developed.  The  accumulation  of  large  internal 
charges  by  the  e-beam  irradiation  at  a  depth  below  the  surface  was  substantiated  by  four  independent 
means;  visible  discharge  channels,  electro-optics,  electro-acoustics,  and  calculation:  Figures  4, 31, 34, 
and  27,28,  respectively.  The  discharge  process  was  studied  in  detail  and  a  three  stage  breakdown 
model  was  formulated.  The  role  of  available  prompt  energy  to  discharge  channel  development  was 
investigated  by  means  of  the  unique  charge  implantation  method.  This  clearly  demonstrated  that  static 
space  charges  can  be  the  source  of  electrostatic  energy  which  drive  a  discharge.  Different  materials 
were  investigated  to  show  major  parameters  which  were  common  and  those  which  were  different. 
The  development  of  dark,  carbonacious  looking  channels  was  found  to  correspond  to  lower 
channeling  resistance  materials.  The  accousdc  and  optic  space  charge  probing  methods  are  new 
diagnostics  which  were  employed  in  this  study  and  are  an  additional  benefit  of  the  program. 
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